Small Angle Scattering data analysis for dense 
polydisperse systems: the FLAG program 



Flavio Carsughi\ Achille Giacometti^ and Domenico Gazzillo^ 

Facoltd di Agraria, Universita di Ancona and INFM Unitd di Ancona, Via Brecce 

Blanche, 1-60131 Ancona, Italy 
Dipartimento di Scienze Amhientali, Universita di Venezia, and INFM Unitd di 

Venezia, S.Marta 2131, 1-30123, Venezia, Italy 
^ Dipartimento di Chimica Fisica, Universitd di Venezia, and INFM Unitd di 
Venezia, S.Marta 2131, 1-30123, Venezia, Italy 



Abstract 

FLAG is a program to calculate the small-angle neutron scattering intensity of 
highly packed polydisperse systems of neutral or charged hard spheres within the 
Percus-Yevick and the Mean Spherical Approximation closures, respectively. The 
polydisperse system is defined by a size distribution function and the macro-particles 
have hard sphere radii which may differ from the size of their scattering cores. With 
FLAG, one can either simulate scattering intensities or fit experimental small angle 
neutron scattering data. In output scattering intensities, structure factors and pair 
correlation functions are provided. Smearing effects due to instrumental resolution, 
vertical slit, primary beam width and multiple scattering effects are also included on 
the basis of the existing theories. Possible form factors are those of filled or two-shell 
spheres. 



PROGRAM SUMMARY 

Title of program: FLAG 
Catalogue identifier: 

Program obtainable from: GPG Program Library, Queen's University of Bel- 
fast, N. Ireland 

Computer: Digital Workstation AU 433 (128 Mb RAM), Pentium I MMX 200 
MHz (64 Mb RAM), Macintosh Powerbook G3 400 MHz (192 Mb RAM) 
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Operating systems: Digital UNIX 4.0E, Windows NT service pack 4, Mac OS 
9.0.2 

Programming language: Fortran 77 

Memory required to execute with typical data: 640 kwords 
No. of bits in a word: 32 

No. of lines in distributed program, included test data, etc.: 6873; of which 
1304 are routines of the Harwell Subroutine Library (HSL) library and their 
use must be acknowledged in any paper publishing results obtained by FLAG. 
The entire code must be linked with the International Mathematical Statistical 
Libraries (IMSL) library. An example of the input file flac. dat is also available 
(94 lines) together with 11 files *.com where common definitions arc found (331 
lines) and a file {flac. par) containing the parameters used for dimensioning the 
arrays (7 lines). 

Additional keywords: Small-Angle Scattering, Hard Spheres, Polydispersity, 
Percus-Yevick, Mean Spherical Approximation. 

Nature of physical problem: The problem is the calculation of the Scattering 
Cross Section in Small- Angle Scattering of polydisperse neutral and charged 
hard spheres. Both dense or dilute systems are considered. 

Method of solution: The algorithms implemented here are obtained by solv- 
ing the Orstein-Zernike integral equations within the Percus-Yevick or mean 
spherical approximation closures for neutral or charged hard spheres, respec- 
tively [1,2]. 

Restriction on the complexity of the problem: Only hard sphere (neutral or 
charged) potentials are used. 

Typical running time: A test run considers 200 points for defining the size dis- 
tribution function and calculates the Scattering Cross Section over 2^^ points, 
without any smearing effects. For neutral hard spheres, on a DIGITAL Work- 
station AU 433 this test takes 7.7 s, while on the Macintosh and on the PC 
26.1 and 49.6 s, respectively. On including instrumental smearing considering 
only one experimental configuration (by convoluting over 1 1 points) the CPU 
time on the DIGITAL Workstation AU 433 increases up to 46.8 s; if multiple 
scattering correction is also taken into account (by using 2'^ points for the 2D 
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Fourier transform), the CPU time is 47.3 s. Furthermore, the addition of a 
vertical sht (sampled by 2^ points) requires 270.1 s of CPU time. For charged 
hard spheres, the calculation without any correction takes 12.1 s of CPU time. 
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1 Small Angle Scattering 

Small Angle Scattering (SAS) is sensitive to the presence of structural and 
chemical inhomogeneities inside the systems, which behave as scattering en- 
tities. Here and in the program, we shall refer to these entities as scattering 
particles whenever they do not bear a charge and macroions otherwise. In the 
latter case a set of smaller scattering cores {counterions with opposite charge) 
are clearly present to satisfy electroneutrality. 

Polydispersity means that the particles are not all identical but their size, 
charge or other properties exhibit a large variety of values (a phenomenon 
commonly found, for instance, in colloidal and micellar suspensions). Hence 
in presence of polydispersity SAS intensity may greatly differ from that of 
monodisperse systems. The aim of the present paper is to present a code 
(FLAC) which can be used in the analysis of experimental data of small-angle 
neutron scattering from polydisperse systems. The FLAC program provides 
the scattering functions from polydisperse fluids of neutral [1] or charged hard 
spheres [2], based on the analytical solutions of the Ornstein-Zernike (OZ) 
integral equations within the Percus-Yevick (PY) or the mean spherical sp- 
proximation (MSA) closures, respectively [3]. For the sake of clarity and in 
order to be self-contained, general formulas appearing in the SAS theories will 
now be briefly recalled. 

1.1 General Equations for polydisperse systems 
The most general expression of the SAS is [4] 
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Here ^(q) is the macroscopic differential cofierent Scattering Cross Section 
(SCS) as a function of tlie exclianged wave vector q, whose magnitude is 
defined by 

47r 

g = — sin ^ (2) 



where A is the neutron wavelength and 26 the full scattering angle. In Eq. 1, 
the integral is extended over the total sample volume V, r is the position 
vector and p(r) the local scattering length density. For X-rays, the scattering 
length density is replaced by the electron density, leaving the general formalism 
unchanged. The angular brackets (• • ■) represent an ensemble average over all 
possible positions and orientations of the particles in the system. 

In the absence of long range order, the scattering length density can be thought 
as having a uniform value po on which fluctuations Sp{r) are superimposed. 
The distribution of scattering material inside a particle of species i is then 
defined by its form factor 

B,{q) = jfdr 6p{r) e'^'^ (3) 



where 

fi^Jdr 5p(r) (4) 

Vi 



is the scattering amplitude at zero-angle, and Vi is the volume of the scatter- 
ing core (i.e the subrcgion of a macroparticle which contains the scattering 
material: nuclei, etc.). Often it proves convenient to combine Eq. 3 and 4 as 

Let us assume that the polydisperse system can be regarded as a mixture with 
p components. By introducing an effective form factor P(q) as 

P(q) = X:n,/^|5,(q)r, (5) 

i=l 



one can also define an effective structure factor S{q) of the mixture by rewrit- 
ing Eq. 1 in the form 

^(q)=P(q)5(q). (6) 
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Here, -P(q) accounts for the distribution of scattering matter inside the parti- 
cles, while S'(q) is related to the particle-particle interactions and the resultant 
equilibrium structure of the system, and it is defined through Eq. (6). 



1.2 Isotropy 



If the particle-particle interactions are spherically symmetric (homogeneous 
and isotropic fluids), then the SCS can be cast into the Fournet-Vrij's form 

[1]: 



1=1 



\{BM)). 



(7) 



where rii is the number density of species i and the brackets {■ ■ ■)u) repre- 
sent angular averages (this formula being still valid in case of non-spherically 
symmetric particles). The Sij{q) are the Ashcroft-Langreth partial structure 
factors [5], defined by 

Sij{q) = 5ij + ,/n~r)~hij{q), (8) 



where 5ij is the Kronecker delta and hij{q) is the three-dimensional Fourier 
transform of the total correlation function hij{r), which can be obtained by 
solving the OZ integral equations of the liquid state [3] . 

The first term on the rhs of Eq. 7 is called the "incoherent" part, and the 
remainder the "coherent" one. The incoherent part vanishes if the distribution 
of scattering matter is spherosymmetric and in this case 

(A(q))^= {B;{q))^^B,{q) (9) 



1.3 Two-Phase Model 



Fortunately, for most of the investigated systems the particles are formed by 
a single type of uniform scattering material (first phase = solute) and are 
embedded in a uniform medium (second phase = solvent for a solution, or 
matrix for an alloy). In this case, Sp{r) in Eq. 3 can be replaced by Ap = Pp—po, 
where Pp and Po are the scattering length density of the scattering particles 
and the embedding medium, respectively. 
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Thus the zero-angle scattering amphtude (Eq.4) is given by 

fi = Ap V, 



(10) 



and the form factor Bi{q) (Eq. 3) becomes 



^i(q) = ^ I dr e'^'' (11) 



Vi 



1.4 Spherical scattering cores 



In this case, Eq. 7 takes the simple form 

— (g) = (Ap)' Y: Vi Vj BM) Bj{q) (12) 



Then the effective structure factor can be expressed in terms of partial struc- 
ture factors as 

c-/ N _ EL■=l^/n^ Fi{q) Fj{q) Sijjq) 



The Bi{q) can be calculated analytically with the result 

Bt{q,R^)= 3^^^, (14) 
qUi 



where Ri is the scattering core radius of a particle of species i, and ji{x) = 
(sin a; — a; cos x) /x^ is the first-order spherical Bessel function. In general, Ri < 
i?f with Rf^ being the hard sphere radius of the i-th particle. Clearly, 
Fi{q,Ri)^ApViBt{q,Ri). 



1.5 Hollow spheres 



In the case of micelles or in the presence of a depletion zone around a pre- 
cipitate, the scattering particles can be thought as a two-shell spheres, char- 
acterized by an inner core and an outer shell with different scattering length 
(or electron) densities. Even in this case the form factor Fj(g) has a relatively 
simple form 

Fi{q, Ra, R,2) = (pi - P2) Va B'^{q, R,,) + (p2 - Po) B-^{q, R,^), (15) 
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where Rn and Ri2 are the radii of the inner and outer shell, respectively, 
Vim = 47r7?f„j/3 (m = 1,2) and Pq, pi and p2 are the scattering length (or 
electron) density of the solvent, of the inner and the outer shell, respectively. 



1.6 Smearing effects 



Experimental data are in general plagued by several effects which might af- 
fect the scattering intensity, as compared to the theoretical one, and hence 
may lead to an incorrect interpretation. Most of them have been well studied 
and the way of controlling them is now well established. The instrumental 
resolution (i.e. the incident radiation wavelength spread, the finite coUimation 
and the detector resolution) together with the radial average of the SAS data 
performed when isotropic data are collected on a 2D detector can be com- 
puted following the recipe given in Rcf. [6], where each effect is described by 
a gaussian contribution, which is also used to model the primary beam for 
the beam width effect [4]. The vertical slit effect is also taken into account 
in the calculation, as described in Ref. [4]; this correction becomes relevant 
for slit systems and double crystal diffractometers. Finally, possible correc- 
tions due to multiple scattering arc also calculated using the MUX routine [7], 
which is based on the theory developed in Ref. [8] . Since SAS experiments are 
performed with more instrumental configurations in order to scan the largest 
possible q range, up to three different setups are considered during smearing 
of the scattering intensities. 



1.7 Backgrounds 



Two different types of background are considered in FLAG. A flat incoherent 
one, like a plateau on which the small angle scattering intensity is superim- 
posed, and the Porod one. The latter is due to the asymptotic scattering of 
large inhomogeneities and can be approximated as [4] 

^yi Porod 

:^ (?) ^ kp q-' (16) 



where kp is the Porod constant. 
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2 Interpeirticle potential 



As often reported in the literature (see e.g. Ref. [9]), the hard sphere radius 
Rf^ of a particle need not coincide with the radius i?j of its scattering core. 
To include this possibility, the ratio s = Rf^/ Ri can be tuned in the program. 



2.1 Neutral hard spheres 

For neutral hard spheres, the interparticle potential Vij{r) is described as 



+00 for r < i?/^-^ 

(17) 

for r > 



where R'^^ = {^f^ + ^f^)- The OZ equations for an hard sphere mixture 
can be solved by using the PY closure. A closed-form for the SAS intensity 
for polydispcrsc hard sphere fluids was obtained by Vrij [1]. Our code for this 
part is constructed along the lines of that reference. 



2.2 Charged hard spheres 



In the case of charged hard spheres, the interparticle potential Vij{r) takes the 
following form 



+00 for r < Rj^,^ 



V^J{r) = { (18) 
"^ZiZj/^er) for r > R^^^ 

where Coulombic (attractive and repulsive) interactions are added to the pre- 
vious neutral case, e being the elementary charge, Zi the valence of the i-th 
species and e the dielectric constant of the solvent. This potential defines the 
Primitive Model of ionic fiuids [3]. The MSA analytical solution of the OZ 
equations for this model allows to write a closed-form for the SAS intensity 
of polydisperse ionic fluids [2]. We assume that polydispersity is associated 
to macroions, whereas counterions are monodisperse. The electroneutrality 
condition, necessary to have a stable system, then reads: 

p 

^ Ui Zi + ric Zc = 0, (19) 

i=l 
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where rtc and Zc are the number density and charge of the counterions, respec- 
tively. 



3 Program description 

3. 1 General 

The FLAG program is able to calculate the SCS of neutral or charged polydis- 
perse hard spheres. In both cases FLAG also provides the result corresponding 
to non- interacting spheres (i.e. P{q))- The scattering density can correspond 
either to a filled or to a two-shell hollow sphere. These models are suitable for 
a large variety of physical systems, among which colloids and mycelles. 

FLAG can work in two different modes: 1) it can provide the theoretical calcu- 
lation of the scattering function, given a set of q (experimental or simulated) 
values ; 2) it can fit experimental data by using the two models described 
above. 

The first mode is particularly useful to gauge the input parameters as a start- 
ing point for the fit. This is governed by the ifit fiag and it will described in 
detail below. 

The filenames to be used are handled in the subroutine file_name and the 
proper input data are read by using the subroutine read_data. Subsequently, 
the subroutine control checks the consistency of the input values. If required, 
the definition of the smearing effects is worked out before the start of the cal- 
culation of the SGS in the fres_prep subroutine. The VA05AD routine of the 
Harwell Subroutine Library (HSL) [10] needs the definition of the functional 
to be minimized in the calfnf subroutine and the statistical errors are esti- 
mated in the kovafd routine developed at the Institut fiir Festkorperforschung 
of the Forschungszentrum Jiilich, Germany. In input FLAG requires a set of 
physical parameters (average size of the particles, volume fraction, etc) plus 
the parameters of the distribution. In addition to the two parameters of the 
distributions considered here (Gaussian, Schulz, Weibull and log-normal) a 
parameter for the normalization is also required. We remark, however, that 
these three parameters are actually computed in terms of the above physical 
parameters. For each q value, the intensity is calculated by using the subrou- 
tine intensity, where the proper model is selected. At this stage, the Porod 
background is considered. In sequence, the corrections due to instrumental 
smearing effects, vertical slit and multiple scattering are performed. Since a 
flat background does not alter any of the previous corrections, the incoherent 
background is added at the final stage. During the fit, a file called *.int (see 
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below) is output at each iteration, for a clearer control of the convergence of 
the numerical procedure. Upon convergence, FLAG estimates the errors in the 
physical quantities and reorganizes the data for the output files. 

3.2 Input/Output description 

At the initial stage FLAG requires a filename interactively and then it sets 
up input and output filenames by appending the proper extension to this 
filename. A short description of each extension is reported below along with 
the meaning of all the files. We note that only the *.dat (containing the 
various parameters) is strictly necessary for running FLAG, the experimental 
data in *.exl, *.ex2 and *.ex3 being necessary only in the case of fit. In 
this case the three files refer to experimental data obtained at three different 
experimental configurations that are are then assembled together into a single 
array. All the other output files are created by FLAG. 

The quantities contained in each file are organized in columns whose labels 
are reported in brackets (used here for clarity and not appearing in the file) 
and have their meaning written underneath. 

Input files: 

dat: general input data 

exl, ex2, ex3: experimental data (Q 1 E) of data collected in different instrumental con- 
figurations. This file is required only in case of fit. 
Q: Scattering vector q (A~^) 
1: Experimental intensity (cm^^) 
E: Experimental uncertainty (cm~^) 

Output files: 
res: General results 

int: Galculated quantities (Q E D G B A) 

A: Theoretical model + Porod background + instrumental resolution (cm~^) 

B: A + slit effect (cm^^) 

G: B + incoherent background (cm~^) 

D: G + multiple scattering (cm~^) 

E: D + random error (simulation) / experimental data (fit) (cm~^) 
Q: Scattering vector q (A~^) 

teo: Galculated quantities (Q D G B A) 
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A: Effective structure factor 

B: Effective form factor (cm^^) 

C: Tlieoretical model (cm~^) 

D: C + Pored background (cm~^) 

Q: Scattering vector q (A~^) 

gen: Calculated quantities on different q range than the experimental one {Qcaic 
C B A). This file is generated only in case of fit. 
A: Effective structure factor 
B: Effective form factor (cm~^) 
C: Theoretical model (cm~^) 
Q: Scattering vector qcaic (A~^) 

fon: Backgrounds (Q P F) 

F: Flat incoherent background (cm~^) 
P: Porod background (cm~^) 
Q: Scattering vector q (A~^) 

nr: Scattering particles (or scattering cores) and hard sphere size distributions 
(R A B D M N) 

R: Scattering particle radius (A) 

A: Scattering particle size distribution (cm~^) 

B: Error of A (cm^'^) 

D: Hard sphere diameter (A) 

M: Hard sphere size distribution (cm~^) 

N: Error of M (cm~^) 

gr: Pair correlation function (R G) 

R: Distance from the center of a particle (A) 
G: Pair correlation function 

We note that at the beginning of the calculation FLAG transforms all dimen- 
sional quantities in suitable powers of Aso that all calculations are performed 
in such units. In the final output all quantities are then restored to their 
original more convenient units. 

3.3 Detailed description 

All relevant parameters are to be specified in *.dat file. First, the spatial range 
{r-min and r-max) of the size distribution function associated to polydisper- 
sity and its number of points (num-pointsJnt) are required. For simplicity, 
the number of input points is given in powers of 2 throughout the code. Hence, 

for instance, 2'^'"'^-^'^^^^'^-^'^^ is the actual number of points of the distribu- 
tion function. The mesh size in reciprocal space Aq is computed from the 
range extension (q^min and q^max) and the corresponding number of points 
(num^points) . 
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The instrumental resolution (for wavelength band and finite collimation) , de- 
tector, radial average and beam width effects are controlled by separate flags: 
iresl, ires2, iresS, ires4 and iresS, respectively. The number of points for the 
convolution (nfold), the radii of source (rl) and the sample (r2) slits, source- 
sample (rll) and sample- detector (rl2) distances, the wavelength (Ham) and 
its band (dlam), the spatial resolution of the detector (rdet), the sampling used 
for radial average (rrav), the FWHM of the primary beam sbeam and the to- 
tal number of experimental configurations nconf have to be provided and they 
will be used only if these corrections are considered in the calculation. Note 
that rll and rl2 are arrays because of the different experimental configurations 
(maximum 3) considered in this analysis. For the vertical slit, the flag (isl), 
the number of points (nsl) used and the maximum q value (tmax) are consid- 
ered as input parameters. The correction for multiple scattering is controlled 
by a flag (ims), and it requires the sample thickness (thick) and the number 
of points to be used for the 2D Fourier transform (nft); for this calculation, 
FLAG is exploiting the IMSL routines [11]. 

The data for the fitting routine (VA05AD) of the HSL are also input parame- 
ters: the step increment (h), the maximum excursion (dmax) and the accuracy 
( acc ) of the numerical calculation can be gauged at will. Moreover, the maxi- 
mum number of iterations (maxfun) as well as the output form (iprint) can be 
also varied ( in this last option, a customized option for a better visualization 
of the evolution of the fit parameters, is also possible). In the general part of 
the input data, one finds the main flag (ifit), which discriminates between the 
simulation of a theoretical curve or the fit of an actual set of experimental data 
(ifit = 1). In the first case the theoretical scattering functions are provided ei- 
ther for a generic range of q (ifit = 2) or for a specific set of experimental q 
values (ifit = 0). This latter case is particularly useful to find a proper starting 
set of parameters for the best fit analysis. 

For both the fit and the simulation of experimental data, one should define 
the first (nl) and the last point (n2) of each experimental scattering intensity 
file to be used and another flag indicating whether the experimental data also 
include the experimental uncertainty (ierr). Note that also nl and n2 are ar- 
rays like rll and rl2. For simulated data, a random error can be introduced in 
the calculated intensity and errjint is the band width where the intensity can 
be randomly shifted using a uniform random number generator controlled by 
the input variable iseed. A weighted fit can be chosen with ipes as well as the 
backgrounds (flat and Porod ones), which are regulated by two flags: ifin and 
ipor, respectively. The flat background can be considered as it stands or it may 
take into account the volume fraction of the scattering particles; this case is 
particularly useful for highly concentrated systems. 

The choice of the size distribution function (idist) allows four different func- 
tions: Gaussian, Schulz, Weibull and log-normal. Galculations for concentrated 
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or dilute systems can be controlled by the inter flag. The form factors of the 
scattering particles can be chosen between the one- and two-shell spheres by iff. 
For charged hard spheres, the counterions charge (zc) and the Bjerrum length 
(lb) defined by 



are to be given. Here kg is the Boltzmann constant, T the absolute temperature, 
e the unit charge and e the dielectric constant. 

Finally a set of n < 13 physical parameters are provided for the calculation, 
each one with its own flag, describing whether its value is kept pinned to that 
given in the input data file ( corresponding flag=0 ), left free for the optimiza- 
tion (corresponding flag = 1) or fixed to a value given in an eocternal file *.res 
produced as output by FLAG (corresponding flag — 2). These parameters are 
the volume fraction xf, the average radius xc of the scattering particles (in the 
case of a two-shell sphere, this correspond to the outer shell radius) and the 
dispersion xd of the particle size; for monodisperse systems xd must he set to 0. 
A condition we have used in our calculation is that the hard sphere radius must 
be larger or equal to the outer dimension of the scattering particle; xs> 1 then 
represents the ratio between the hard sphere radius and the outer dimension 
of the scattering particles. The Porod constant (xp) and the flat background 
(xfo) are required and might also he optimized. The differences of scattering 
length densities between the outer and inner shell (xh), between the outer shell 
and the solvent for the scattering particles (macroions) (xcon) and between 
the counterions and the solvent (xconl) are also necessary for the calculation. 
The thickness of the outer shell for two-shell spheres is set by the variable xa 
and xz represents the charge associated to the macroions with sizes equal to 
the average of the distribution. Finally, a scaling factor xk is necessary when 
data are not available in absolute units. At the end of the input file, the name 
of the external file where some data have to be used is also requested. All the 
input data are tested for non- appropriate values. 

All adopted distribution functions depend upon three parameters. For three 
distributions (Gaussan, Schulz and log-normal) , they can he analytically com- 
puted from the knowledge of the volume fraction, the mean radius and the 
polydispersity index ( dispersion ) which is a measure of the strength of polydis- 
persity and whose definition can be found in Ref [12]. In the Weibull case, the 
lack of a one-to-one correspondence in the analytical expression enforces a nu- 
merical self- consistent procedure. All the n physical parameters to he optimized 
are constraint to be positive, apart from xs which must be larger than 1 and 
xb which can assume all real values. In the case of charged hard spheres, the 
population of counterions is determined by the electroneutrality condition and 
the value of the screening parameter V is computed iteratively using Eq. 50 of 
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Ref. [2]. 



The coherent transmission factor T^oh is calculated by estimating the theoretical 
scattering cross section over a very wide interval in the q space 




where D is the sample thickness and the term between square brackets repre- 
sents the macroscopic coherent scattering cross section S under the approxi- 
mation of isotropic scattering, i. e. the integral of the differential one over the 
solid angle. This transmission factor is an indicator of the multiple scattering 
probability: the higher Tcoii, the smaller is the multiple scattering probability. 
This effect is incorporated in FLAC, by the use of the MUX routine [7]. The 
output of the routine is adapted to the experimental q grid by using linear 
interpolation, with the exclusion of the last point which is obtained by extrapo- 
lation. The structure factor for a system containing monodisperse counterions 
and macroions can also be computed with the charged monodisperse option. In 
this case the scattering is reduced to the sum of 4 terms, corresponding to a 
2x2 matrix. 

The average pair correlation function g{r) is calculated by using the real space 
mesh size Ar determined directly from the one in reciprocal space by using 
the standard relation Ar = 7t/{N Aq), where N is the number of points 
num_points, and it exploits a standard Fast Fourier Transform (FFT) algo- 
rithm. 

3.4 Compilation 

FLAC must be placed in the same directory of the files containing the parame- 
ters (flac.par) and the common definitions (distribution.com, flag.com, gen- 
eral. com, ggc.com, intensity.com, library.com, names.com, parameters.com, 
points, com, quantities.com and smearing.com) and it must be linked with IMSL 
library [11]. 



4 Experiment 

In order to illustrate the practical use of FLAC, we provide, along with the 
code, a set of so far unpublished experimental data collected, by us, at the 
point geometry KWSI diffractometer at the FZJ in Jiilich, Germany, on silica 
particles in hydrogenated water. Samples have been prepared from commercial 
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Ludox HS30 at the nominal volume concentration of 16.5%. The measurements 
have been performed by using 4-56 A neutrons with a wavelength band of about 
20% and the data have been collected with three different experimental configu- 
rations, i.e. sample-to- detector and collimation distance equal to 14, 4 o^i^d 1.2 
m, in order to investigate the largest possible q range [13]. The samples were 
contained in quartz cell of 2 mm thickness, the multiple scattering being then 
not negligible. Hence, we have performed the data analysis by using the neutral 
one-level hard sphere model with 6 free parameters, i.e. the volume fraction, 
the average radius, the dispersion of the size distribution function, the hard 
sphere factor, the incoherent background and the contrast. The Weibull size 
distribution function has been selected in the present case. We have optimized 
the free parameters by including in the calculation also the smearing due to the 
wavelength band, different finite collimations and multiple scattering. In Fig.l, 
experimental data collected at three different instrumental configurations are 
depicted together with the theoretical scattering cross section, the incoherent 
scattering, the effective form factor and the smeared scattering cross section. 
Error bars are smaller than the size of the symbols and hence are not displayed. 
The agreement between experimental data and calculated cross section is ex- 
cellent and the physical parameters shown in Tab. 1 are well representative of 
the system. The coherent transmission Tcoherent is relatively high, showing that 
the coherent multiple scattering does not play a major role in this sample. A 
nominal value of the scattering contrast of commercial Ludox HS30 has been 
estimated to be about 4-10^^ cm~^, in good agreement with our optimized value. 

Acknowledgments 

The Italian MURST (Ministero delVUniversita, e della Ricerca Scientifica e 
Tecnologica) , the INFM (Istituto Nazionale per la Fisica della Materia) are 
gratefully acknowledged for financial support. We thank Raffaele Della Valle 
for a critical reading of the manuscript. The authors are indebted to the HSL 
management for the free use of the routines, and to the computing group of the 
Institut fur Festkorpeforschung of the Forschungszentrum Jillich, Germany, 
for their assistence during the preparation of this work. The DuPont Company 
is also gratefully acknowledged for providing the Ludox HS30 material. 



References 

[1] A.Vrij, J. Chem. Phys. 69, 1742 (1978); 71, 3267 (1979). 

[2] D. Gazzillo, A. Giacometti and F. Carsughi, J. Chem. Phys. 107, 10141 (1997). 

[3] J. P. Hansen and I. R. Mc Donald, The Theory of Simple Liquids, Academic 
Press, London (1986). 



15 



[4] O. Glatter, in Small Angle X-ray Scattering, Eds. O. Glatter and O. Kratky, 
Academic Press, London (1982). 

[5] N.W.Ashcroft and D.C.Langreth, Phys. Rev. 156, 685 (1967). 

[6] J. S. Pedersen, D. Posselt and K. Mortensen, J. Appl. Cryst. 23, 321 (1990). 

[7] M. Monkenbusch, J. Appl. Cryst. 24, 955 (1991). 

[8] J. Schelten and W. Schmatz, J. Appl. Cryst. 13, 385 (1980). 

[9] J. S. Pedersen, Phys. Rev. B47, 657 (1993). 

[10] Harwell Subroutine Library, Release 11, AEA Technology (1993). 

[11] IMSL, Fortran routines, Visual Numeric (1994). 

[12] B. D'Aguanno and R. Klein, J. Chem. Soc. Faraday Trans. 87, 379 (1991). 
[13] F. Carsughi, A. Giacometti and D. Gazzillo unpublished. 



16 



FIGURE CAPTIONS 



Fig.l: Best fit of commercial Ludox HS30 with neutral one-level hard sphere. 
The experimental data collected at I4, 4 1.2 m are shown together with the 
theoretical intensity, the effective form factor, the incoherent background and 
the best fit obtained by smearing for the wavelength band, the finite collimation 
at different instrumental configurations and the multiple scattering. 
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20.3± 0.7 
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81. 1± 0.7 
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HS 




1.22± 0.01 


IB 
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0.954± 0.003 
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(cm~^) 


3.43±0.05-10i° 
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(cm"'') 


8.1± 0.3-10^*^ 


Sp 


(cm"-'-) 


6.9± 0.3-105 


'^coherent 


(%) 


92.9 



Table 1 

Physical parameters obtained by the best fit with a Weibull size distribution function 
of the commercial Ludox HS30. The parameters optimized by the best fit procedures 
are: Cy the volume fraction, R the average radius, ^ the dispersion of the size 
distribution function, HS the ratio of the radii of the hard sphere and the scattering 
particle, IB the incoherent background, Ap the scattering contrast. Moreover, Np 
the scattering particle density and Sp the specific surface are then calculated. The 
coherent transmission Tcoherent aiso estimated. 
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